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Introduction. Graphical representations have provided a powerful tool for both the analyt- 
ical and numerical study of statistical mechanics systems. However, the full power of these 
methods can only be realized when the phases and the phase transitions in the original system 
are easily identifiable in the graphical representation. Usually this involves the percolation 
properties of the graphical representation. Consequently, the subject is still something of an 
"art" ; various systems are treated on a case-by-case basis. Despite the recent progress on sys- 
tems other than Potts models (c.f. [CM II] and references therein) an apparent requirement 
is a high degree of internal symmetry for the individual degrees of freedom. 

In this paper, we will treat the simplest possible system with no internal symmetry, namely 
the ferromagnetic Ising model in an external field. The representation is actually that of a 
duplicated (or "replicated") system and is closely related to the so called red/blue represen- 
tation of the Ashkin-Teller model that was introduced in [CM I]. Before launching into a 
full description of this representation - the subject of the next section - let us illustrate the 
problems that beset the standard approaches. 

Consider the usual Ising Hamiltonian: 

h= — JijQ-jQ-j - y^^hjQ-j (i) 

where the bonds ((i,j)) and sites (i) are part of some regular graph, the Jij are positive and 
— oo < hi < +oo. To simplify a number of technical points, it will be assumed that the model 
is short ranged, that is to say, the number of bonds emanating from any site is uniformly 
bounded. 

In case hi = 0, there is the standard random cluster representation [FK] for bond configu- 
rations uj: In finite volume let 

W(oj) = B 3 (u)2 c ^ (2) 

where Bj is the Bernoulli factor of 1 — e~ 2,3JiJ for each occupied bond of uj and e~ 2,3JiJ 
for each vacant bond and c(ui) is the number of components. (Here we ignore modifications 
due to boundary conditions.) Since, for each u, these quantities (weights) are non-negative 
they define a probability measure on the collection of all bond configurations. As such, these 
measures admit infinite volume limits. It is remarked that for regular infinite volume systems, 
e.g. nearest neighbor models on Z d , the percolation of these bonds corresponds to positive 
spontaneous magnetization in the spin-system [CK, ACCN]. 

For hi not identically zero, there are two obvious approaches both of which lead to un- 
satisfactory results: First we may include the magnetic fields directly in the formula for the 
weights. The result is 

C ( W ) 

W(u) = Bj (uj) 11 2 cosh[h(K Q )] (3) 

a=l 

where the product runs over the clusters Ki, . . .K c ^ and h(K a ) = (3J2 ieKa hi. In other 
words each cluster is given a weight in accord with the field energy obtained if all spins of the 
cluster are of the same sign. The second approach - which is most useful if all the (hi) are of 
the same sign - is the introduction of ghost bonds. In general, the ghost system is defined as 
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follows: clusters are defined and preliminarily weighted as in the zero field case. In addition, 
for every site i such that hi > 0, there may be an "up" ghost bond which carries the weight 
I — e~ 2 ^ hi or no such bond with weight e~ 2 @ hi . Similarly, if hj < 0, the site j may get a 
"down" ghost bond with weight 1 — e _2/3 '^ J ', etc. However, the system is constrained so that 
no cluster may have ghost bonds of both types. Thus, if To denotes the extended system of 
ordinary and ghost bonds the weight for To is given by 

C( W ) 

W(TJ) = <f>p)Bj, h p) J] v a - (4) 

a=l 

here $(uJ) is the function that is one if the above mentioned constraint is satisfied and zero 
otherwise, v a is 2 if there are no ghost bonds in the a th cluster and 1 otherwise and Bj^(u) is 
the Bernoulli weight for the configuration uJ: 1 — e 2f3JiJ (e 2 ^ Ji - J ') if the bond is occupied 
(unoccupied) and l — e 2f3 \ hi \ (e 2 ^^) if the ghost bond at i is occupied (unoccupied). It is not 
hard to show that the weights in Eqn. (3) are obtained from those in Eqn. (4) by summing 
over the ghost bonds. 

The difficulty with either of these representations is manifestly apparent. Consider the 
case hi = h > on the usual <i-dimensional lattice with d > 2. By the Lee- Yang theorem, 
there is no phase transition as temperature is varied. Nevertheless, it is clear that in the 
graphical representation, there is a percolation transition at some finite Pf(h) < (3 C . Thus, 
these representations are susceptible to false percolation transitions. In less trivial situations 
such as random field problems, it can be argued that such false percolations will also occur. 
Cluster algorithms based on this representation are ineffective because of the weighting of 
cluster flips and the problem of false percolation: Notwithstanding the fact that large-scale 
clusters may be of no particular interest, they tend to become dynamically frozen. For 
example, in the second representation, no cluster can flip unless it is devoid of ghost bonds 
which occurs very rarely for large clusters. These effects lead to severe slowing down both for 
critical and non-critical values of the parameters. 

Attempts have been made to circumvent these difficulties by restricting the size of clusters 
(e.g. [NB] ) - with some degree of success - but in principal, there can be no percolative signal 
of the phase transition in any such representation. 

It is therefore of interest to produce representations for spin-systems in arbitrary external 
fields; the subject of this work. In this regard, we have recently learned of a related approach 
to these sets of problems that is applicable to systems with periodically varying fields. These 
methods are described in [HBb] and [HBc] and were applied to the hard-core lattice gas (a 
limiting case of the staggered field problem) [HBa], [HBb] and related systems [HBc]. We shall 
describe the similarities (and differences) between this method and the current method in the 
conclusions section. At present, these approaches are, by and large, restricted to Ising-type 
systems. Nevertheless, this covers a number of problems that are of current interest such as 
the random-field Ising model. In the context of the staggered-field Ising system, an algorithm 
based on this representation has already been tested [RMC] with quite satisfactory results. 

Red/Blue Representation in Zero Field. Let us start with the zero field case. Much of 
this material appears in [CM I] (Appendix B3) in the context of the Ashkin-Teller model. 
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Consider the duplicated system so that at each site i, we have Oi = ±1, Tj = ±1. We will 
envision the duplicated problem as a 4-state spin model, in particular, the 4-state clock 
model. Let us use Si to denote the spin variable, regarded as one of four points on the circle 
and s as notation for a spin configuration. If is a bond of the lattice, there are three 
possible energy states for the spins Sj and Sj, the highest of which is set to zero without loss of 
generality. Thus we define • = 4Jij and = 2Ji,j- The Hamiltonian for the duplicated 
system is given by 

H(s) = - Y^KjSs i=ai + SljS ai=ai+1 + £&.<=.,-!] (5) 

(i,3) 

where the 5's are one if the indicated condition is satisfied and vanish otherwise and we have 
used additive notation for the group structure of the spin-states. Explicitly, s — > s + 1 means 

(++) - (+") - (— ) - ("+) - 

For each we may expand the statistical weight associated with H as: 

e < ) «. i =. j ^«. i =. j +i^i. i =. ) -' = 1 + B itj S ai=a . + Ri,jS Si=Sj+1 + s , , (6) 

with Bij = e@ £i 'i — 1 and Rij = e^ £i ^ — 1. Note that no terms involving products of 5's 
appear since the conditions are exclusive. Thus, there are four possibilities for each A 
"vacant" bond corresponding to the 1, a "blue" bond corresponding to the term Bij5 Si=Sj 
and one of two types of "red" bonds which forces Sj = Sj ± 1. 

A graphical configuration r\ is obtained by a selection of one of the above four terms for 
each bond of the lattice. The weight of r] is given by 

H!i)=En%)fc(p) ( ? ) 

S. <M> 

where (77) = 1, Bij or Rij according to whether a vacant, blue or red bond has been 

selected and 5{ij)(rj,s) is the corresponding delta factor (or 1 for the vacant bonds). 

These weights are not particularly difficult to evaluate explicitly nor are the results partic- 
ularly illuminating. Of greater interest is the weight associated with all configurations r/ with 
the same colored bonds but different types of red bonds. We will denote these equivalence 
classes by an unadorned rj and define V(rj) = Yl v eri—^' ^ n ex P nc ^ formula for V(rj) is also 
possible which we will give after the relevant notation is developed. 

First it is noted that not all 77's have non-zero weight. In particular, it is necessary and 
sufficient that each closed loop of non-vacant (hereafter "occupied") bonds of 77 must consist 
of an even number of red bonds. Let D(r/) be the indicator that the above condition is 
satisfied. Next, let F{rf) denote the "bond fugacity" prefactor (which is the same for all 77's 
in the equivalence class) 

f(i)=n%(i)- ( g ) 

(id) 

1 We remark that the Hamiltonian in Eqn. (5) differs from (two copies of) that in Eqn. (5)) - with hi = 
- by an overall constant 
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Finally, the geometric factors: It is observed that there are two relevant notions of connect- 
edness in a configuration rj, namely connected by blue bonds or connected by red and blue 
bonds, which we will call a grey connection. Let Cb{rj) denote the number of blue components 
in the configuration rj and K(rj) the number of components by the less stringent (grey) def- 
inition. Thus K(rj) < Cb(rj); for example, if each bond is occupied and red, K(rj) = 1 while 
Cb{rf) is the total number of sites. 

With the above in hand, the result is 

V{rj) = F(rj)D(ri)2 K ^2 Cb ^ . (9) 

We remark that we have not yet discussed boundary conditions; the above is strictly true only 
for free boundary conditions on an arbitrary (but finite) graph. We denote the associated 
measure by v(-) with various super- and sub-scripted adornments to be appended as needed 
to indicate the lattice, the boundary conditions and the fields. 

We introduce another measure, called the Edwards-Sokal measure, for the representation. 
This is a joint measure on bond and spin configurations that assigns an equal weight to each 
spin configuration that is "consistent" with the given bond configuration: 

V ES (s, V ) = F(r ] )A(s, V ). (10) 

In the above, A(s, rj) vanishes unless (a) the spin values are constant in each blue component 
of r\ and (b) the spin values Si and Sj across each red bond of rj differ by one: Sj = Sj ± 1. 
If (a) and (b) are both satisfied then A(s, rj) = 1. The Edwards-Sokal measures defined by 
these weights will be denoted by v ES (— ). It is not difficult to see that 

J2v ES (s, V ) = e-P H & (11a) 
v 

and 

J2V ES (S,V) = V( V ). (lib) 

s 

Thus the bond marginal of v ES (— ) is the red/blue random cluster measure and the spin 
marginal is the Gibbs measure for the system. The Edwards-Sokal measure facilitates an 
interpretation of the bond configuration rj. In particular, connected clusters of blue sites 
represent spins that are all in the same state, separate clusters not connected by any bonds 
are uncorrelated from one another and sites connected by an even (odd) number of red bonds 
represent spins that agree mod 2 (differ by ±1). 

Boundary conditions are readily accounted for in the red/blue random cluster measures. 
If A is the graph and <9A the boundary sites, a single spin specification on <9A divides <9A into 
(at most) four components. These components are treated as a single site in the counting 
of K(r]) and Cb(rj). Furthermore, there are restrictions on the configuration rj, e.g. no two 
of the boundary components can be connected by blue bonds. Of particular interest are the 
wired boundary conditions where all the boundary sites are treated as a single component. 
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In the absence of fields, the boundary type may be identified as any of the four states or 
superpositions thereof. 

We now describe a cluster algorithm based on the red/blue representation. Starting from 
a spin configuration s, each is checked for the conditions Sj = Sj or Sj = Sj ± 1. If 

Si = sj, a blue bond is placed on with probability bij = Bij/[1 + Bij] = 1 — e~ 4:f3JiJ or 
is declared vacant with probability 1 — bij. If Sj = Sj ± 1, a red bond is placed on with 
probability r^j = Rij/[1 + Rij] = l~e~ 2f3JlJ or the bond is declared vacant with probability 
1 — rij. On all bonds for which Sj and Sj are completely opposed (sj = Sj + 2, i.e. Si = —Sj) 
the bond is always vacant. All bond moves are performed independently. This creates the 
configuration rj. Given the configuration rj, the updated spin configuration is created by two 
types of cluster moves. First, the spins in each isolated grey cluster can be coherently rotated 
by 0, ±1, or 2 with equal probability. Second, the spin type of all spins in a blue cluster 
(including single spins with no blue bonds) can be reversed (s — > — s) with probability 1/2. 
We will refer to these as grey moves and blue moves respectively. The cluster moves are 
all performed independently. This algorithm is of the Swendsen-Wang [SW] type and it is 
straightforward to show, following [ES] , that it simulates the Edward-Sokal measure of Eqn. 
(10). For more details (and more formal descriptions) concerning general algorithms of this 
sort, the reader should consult [CM I], [CM II] 

Remark. The strength of this representation is that percolation (essentially of blue bonds) 
can be directly related to the non-vanishing of the spontaneous magnetization. This will be 
discussed in full detail in the next section. One principal difficulty is the presence of large 
grey clusters which, under a variety of conditions, will hinder rapid equilibration. Indeed, 
at moderate temperatures above the magnetic ordering temperature, there is, generically, 
percolation of grey bonds. (For example, this can be established by rigorous bounds for the 2d 
square lattice.) Furthermore at low temperatures, the rapid "tunneling" events characteristic 
of better representations are inhibited. For example consider a situation where j3 >> 1 and 
suppose that all the spins on the left side of a box are in the state ++ while on the right 
side of the box all spins are H — . Here both the left and right sides will be tied up in a single 
cluster. If we use two independent SW algorithms (for the Ising case, S° = 2S 1 or, in the 
general Ashkin-Teller case the bilayer representation) the two clusters are disjoint and the 
interface will disappear in a few Monte Carlo steps. However with the red/blue representation, 
there would nearly always be red bonds connecting the left and right regions which prevent 
full alignment of the box; in this case, an extremely long time is required for equilibration. 
This represents the principal reason that this approach was abandoned as a tool to study the 
Ashkin-Teller model. However, as we shall see, in the presence of external fields, the red/blue 
representation has significant advantages over the alternatives. 

Red/Blue Representation with Fields. Before we actually introduce the fields, let us 
present one more refinement which, in the absence of fields is somewhat redundant. Consider 
a configuration rj. Recall that each blue cluster of r\ represents spins that are all in the same 
state. Now a blue cluster may be part of some larger grey cluster but the dynamics described 

in the previous section allows the spins of each blue cluster to flip (++ <-> , H — <-> — V) 

leaving the remaining spins of the grey cluster invariant. But then a further classification 
is allowed for the blue clusters namely whether the cluster is in the (++, ) or ( — h,H — ) 
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states. In the former case, we will say that the cluster is dark-blue and in the latter, light-blue. 
As we shall see, percolation of the light-blue bonds represents the onset of magnetic ordering. 

We will denote these red, dark-blue and light-blue (RDLB) configurations by S's. It is 
seen that the weight of a configuration S is given by 



where F(E) has its previous meaning, Cb(S) is the number of blue clusters - both dark and 
light - and D(S) indicates that within each grey cluster, all paths between pairs of light- 
blue clusters consist of an even number of red bonds, similarly for pairs of dark-blue clusters 
while any path between a dark-blue and a light-blue cluster consists of an odd number of 
red bonds. We note that the "grey moves" of the algorithm described in the previous section 
(where the spins of the grey clusters were coherently rotated) consists, with probability 1/2, 
of exchanging dark-blue for light-blue within each grey cluster. As such it is easy to recover 
the weights in Eqn. (9) from those of Eqn. (12) - for each grey cluster, one sums over both 
these possibilities. We denote the measure defined by the weights in Eqn. (12) by V(— ) and 
the associated Edwards-Sokal measure by Y ES (— ). 

Finally we remark that the division of blue into dark and light degrees of freedom creates 
more candidates for the boundary conditions. In particular we must now separately consider 
light-wired boundary conditions corresponding, in the spin system, to H — or — h boundary 

conditions or dark-wired corresponding to ++ or boundary conditions. We reiterate that 

in the absence of external fields, these are all equivalent: the graphical problem is invariant 
under the exchange of dark and light labels. 

We are finally ready to introduce the fields. Consider the Ising Hamiltonian as given in 
Eqn. (1). Upon duplication we have 



The result is clear: the fields do not act on the spins Si that are in the states (H — ) or ( — h) 
and hence on the light-blue clusters. In the construction of the graphical representation, we 
may either go the route of ghost bonds (as in Eqn. (4)) or field energies (as in Eqn. (5)). The 
latter choice is somewhat easier and gives us the weights 



V(H) = F(E)3(E)2 Cb ^ 



(12) 



H = H(s) = H(a) + H(t) = - ^ J l , J [a l a 3 + t iTj ] - ^ h^a, + n}. 



(13) 



C* db (H) 

V(H) = F(S)©(S)2 C ^ S) J] 2cosh[2h(c*)] 



(14) 



where is the a th dark-blue cluster of H and h(c^, 6 ) is defined by 




(15) 



We again emphasize that these weights are non-negative and therefore define a probability 
measure on the space of all RDLB configurations. As above the measures will be denoted by a 
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V; in particular for the Hamiltonian H at inverse temperature /?, in volume A with boundary 
condition * we will denote these measures by Ya*-/3,h- 

A cluster algorithm based on the RDLB representation is similar to the algorithm described 
in the previous section. The bond moves are identical. There are now three kinds of spin 
moves; grey, dark-blue and light-blue. The dark-blue cluster flips occur with probabilities 
weighted by the exponent of their field energies: For the cluster c^ 6 the ratio of forward to 
backward rates is e ±2h ( c « ) with the sign depending on the initial and final spin state. The 
grey cluster move exchanges light and dark blue cluster with a probability that depends on 
the fields. Although the grey and dark-blue moves depend on the fields in a complicated 
way, the light-blue move is independent of the field. The crucial point is that since the field 
energy is always zero for the light-blue clusters, these are "freely" flipped with probability 
1/2. As we demonstrate in the next section, for all ferromagnetic Ising problems, the light- 
blue infinite cluster density can be identified with the magnetic order parameter. Thus, in 
this representation, the degrees of freedom contributing to the long range order are decoupled 
from the fields. Nevertheless, the algorithm may experience slowing due to the presence of 
large grey or dark-blue clusters. In the final section, we will discuss some strategies that can 
be used to overcome these shortcomings. 

The order parameter. In the following we present some elementary results for Ising sys- 
tems. Although most of these are proved explicitly or implicitly in the standard treatments 
(e.g. [G]) for completeness we will provide a proof sketch. Below we will assume that A is a 
finite graph that is a subset of some infinite L and the Ising Hamiltonian, 7i is as given in 
Eqn. (1). We will consider sequences (A&) which "tend to L" meaning A^ + i D A& and each 
i G L is, eventually, in some A&. This is denoted by A / L. Thermal states on A at inverse 
temperature j3 and boundary condition # on L \ A will be denoted by (— )^ H . The special 
(and most important) cases are # = +(# = —) corresponding to each spin of L \ A set to + 
(-)• 

Proposition 1. Let i G L and for each A with i G A let 
and 

Then the limits 

mf = lim mf(A) 

and 



m- = lim m- (A) 

1 A/X 1 

exist for any sequence (A&) as described above and are independent of the sequence. Let 
Hi((3) = \\rrii(P) — m ^{P)] denote the "magnetization excess". Then the necessary and 
sufficient condition for Gibbsian uniqueness is that for each i, Hi((3) = 0. 

Proof. We use the fact that the finite volume measures (and hence the infinite limits thereof) 
have the strong FKG property for the ordering induced by (<7j = +) > (c^ = — ). This means 
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that the + boundary conditions on A are "higher" than any other boundary conditions. In 
particular, if A' D A the restriction of (-)X'-p to A is dominated by (-)X-p n- This ensures 
the existence of an unambiguous limiting plus state (independent of the sequence (A&)) which 
will be denoted by (— )p n as well as the existence of the mf. Similar considerations provide 

the infinite volume state (— )p- H and the m~ = (o~i)p n - If, for any i, 2^ = mf — m~ is not 
zero then obviously there is more than one state. On the other hand if for every i, fii = 
then, by Strassen's theorem, [S] we have (— )~p n = (-)p'n- smce an y limiting state lies 
above (— )p~ n and below (—)p n (in the sense of FKG) it follows that the limiting state is 
unique. □ 

Remark. Under the usual physical assumptions that L = Z d or some other regular <i-dimensional 
lattice and that the hi and Jij are periodic or i.i.d. random variables then 




-m-(A)} 



exists (almost surely in the random cases) for all A /* L in the sense of Van Hove. The 
quantity /u(/3) may be identified with a thermodynamic derivative. Indeed if f(b) is the free 
energy in the presence of the applied field b (that is to say the term — J2i b&i is added to H) 
then = -§/*[/'(& = 0+) - f(b = 0")]. 
Our main result: 

Theorem 2. Let A denote a finite graph (which may be regarded as a subset of some infinite 
h) and consider the light-blue wired measure Ya™-,p,h(—) on RDLB configurations as described 
in Eqn. (13)-(15) with wired boundary conditions. Let i 6 A and let {i <-> <9A} denote the 
event that i is connected to the boundary by light-blue bonds. Then 

V AW ({z <- dA}) = l[m+(p,A) - m-(p,A)} = ^(/3; A). 

Proof. Let (-)X^p h denote the thermal state for the (four-state spin) model on A with Hamil- 
tonian H as in Eqn. (13) and H — boundary conditions. It is clear that if A is any function 
of the cr's alone then (A)^^ H = (A)~j^.p n (because the a's and the r's are independent) and 
similarly if A' is a function of the r's alone then (A') J ^Jj 3 H = (A')^.p n . Thus 

- n])tp,H = ^K + (A) - m-(A)] = ^(A). (16) 
The stated identity follows easily from the Edwards-Sokal measure. We start with 

- n])Xr^ H = ^ E J w {\[o-i - n}) (17) 

where E^£(— ) denotes expectation with respect to the measure V^£(— ) with light-blue 
wired boundary conditions on dA. Let us consider the three possibilities for the site i: (1) 
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the event {i <-> <9A} occurs, (2) the site % belongs to a dark-blue cluster and (3) the site % 
belongs to a light-blue cluster that is not connected by light-blue bonds to the boundary. 
We may express the right hand side of Eqn. (17) in terms of the three associated conditional 
expectations. 

Whenever the third possibility occurs, regardless of whether % is grey connected to <9A, the 
cluster of i is H — and — h each with probability 1/2. If i is in a dark-blue cluster then cxj = r^. 
Thus, items (2) and (3) contribute nothing. On the other hand, if {i <-> dA} occurs then the 
boundary conditions force Oi = + and n = —. This establishes the stated identity. □ 

Remark. It is observed that this representation may also be used for the calculation of corre- 
lation functions. For example, consider an Ising system in finite volume A with fields (hi) and 
boundary condition *. Duplicating the system with the same boundary conditions it is not 
hard to see that the expectation of (uj — Ti)(o~j — Tj) is exactly twice the truncated spin-spin 
correlation function between the sites i and j. Now if the site i belongs to a dark-blue cluster 
- which is necessitated if it is connected to the boundary - then the contribution is zero 
because o~i = 7$. Similarly if j belongs to a dark-blue cluster. Furthermore, if i and j belong 
to separate light-blue clusters there is cancellation because for fixed Tj and o~j, the quantity 
(o~i — T~i) takes on the value 2 and —2 with equal probability. Finally, if i and j are in the same 
light-blue cluster, the result is four. Thus we are left with another identity: The truncated 
spin-spin correlation function is half the probability that the relevant sites are in the same 
light blue cluster. 

With the exception of a few technical points the picture is complete. Given H — (i.e. 
light-blue wired) boundary conditions, there either is or isn't percolation of light-blue bonds. 
If there is, the stated order parameter is positive and if not, there is no ordering of the 
stated type. In fact, these results apply to other systems of this type. For example, the 
previous sentence holds generally for the Ashkin-Teller model in an external field where 
< E\a < However, specific to the Ising system is the stronger statement that if the 

magnetic excess vanishes (n((3) = 0), there is no long-range order of any type. There is a 
second issue that is not unrelated to the first. We have implicitly "defined" percolation to 
mean percolation in the state that ultimately produces the order parameter. (Of course, in 
simulations these are rarely the boundary conditions that are actually used.) What has not 
been ruled out is the possibility that there is light-blue percolation in some limiting state but 
not in the limiting H — state. Although there can be little doubt that in general this does not 
occur, for Ising systems we have a complete proof. 

We start with a precise definition of (light-blue) percolation: 

Definition. Let A C L be a finite lattice and consider the various measures V A # :/3 (— ) on 
RDLB configurations with boundary conditions For i e A, consider the probability that i 
is connected to dA by light-blue bonds maximized over all possible boundary conditions: 

Pi(P; A) = maxV A#;/ ^({z <-> dA}). (18a) 



We define 



P l ((3;oo) = limP,(/3;A). 



(18b) 
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(It is not hard to show that this limit exists independent of the sequence (A k ).) If Pi(/3; oo) 7^ 
for any i, we say there is percolation and if Pi(f3; 00) = there is no percolation. 

For Ising systems, we have the following characterization: 

Theorem 3. For an Ising system as defined by the Hamiltonian in Eqn. (1) there is a unique 
limiting state if and only if there is no percolation. 

Proof. By Proposition 1, there is uniqueness iff (ii = for all i. Obviously Hi((3) < Pi((3; 00). 
It remains to be established that ^i{f3) = =>- Pi(/3', 00) = 0. Here an argument similar to 
that used in [CM II] Theorem 3.8 may be applied. Let * denote any boundary conditions 
that optimize Va* ; /?({« <-> dA}) which, without loss of generality, is taken to be a single-spin 
specification. There are (at most) four components to OA: <9A ++ (*), dA~~(*) and dA^(*) 
corresponding to the four spin-types. We may write Pi((3; A) = P^ (/?; A) + P i l ~(/3; A) with 

^ + -(/3;A)=V A , ; ^({^M+-(*)}) (19) 

etc. By a^r symmetry, it may be assumed that P^ ((3; A) > P i A). 

We claim that in the light-blue wired (i.e. H — ) boundary conditions, the probability of 
{i <-> <9A+~(*)} is at least as large as it is in the * boundary conditions and hence, after the 
infinite volume limit, Hi((3) > |Pj(/3;oo). Let us consider the situation from the perspective 
of the Edwards-Sokal measure. Decomposing P^ (/?; A) according to spin configurations we 
have 

Va* ;/3 ,h({z <-> M+-(*)}) = X)vf.^fe,T)Vf.^({i - dA+-(*)} I a,r) (20a) 

a_ 

The right hand side is of the form of a thermal expectation of a function of spin configurations, 
namely Vf^ H ({i <- dA++(*)} \ g_, r); 

V a * ;/3 ,h({z - M+-(*)}) = (Vf^ )H ({z - 0A+"(*)} I z,r)>^. A . (20b) 

However, in light of the algorithmic construction of the configurations S given a it is clear that 
this function is increasing in a and decreasing in r: Indeed, this is amounts to a bond-site 
percolation problem and increasing the number of plus cr's and minus r's just increases the 
number of sites. But then, by the FKG property of the Gibbs measures, 

(Vf*V({* ^ <9A+»} I a, r)) A;/3jH < (Vf£ ;/Mr ({i - dA + ~(*)} \ a, t))+^ >h 

= V A^ , ;/3 , ^^ ({^ ^ M+-(*)}) < ^(/3; A) (21) 

from which the stated result follows. □ 

Conclusions. The upshot of the previous analysis is that the crucial degrees of freedom at 
points where there is a magnetic ordering transition are represented by the light-blue clusters. 
For the algorithms described, the spin-transitions within these clusters is uninhibited and this 
should lead to substantial improvement over other cluster methods in many situations. On the 
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other hand, when the system is critical or at a point of phase coexistence, it may be too much 
to hope that there are never any other large-scale fluctuations. In particular, in the absence 
of fields, by light/dark symmetry, the distribution of dark-blue clusters is identical to that 
of the light-blue. In the presence of fields, there might still be large dark-blue clusters and 
these would have a tendency to "freeze" . Furthermore, since we expect grey percolation at all 
points of physical interest, the transitions that exchange dark with light could be disastrously 
slow. 

However even assuming the worst case scenario, namely that the only aspects worth pre- 
serving are the light-blue moves, the algorithm may be supplemented with Monte Carlo steps 
from other algorithms to insure ergodicity. In these circumstances, the whole of the previously 
described algorithm need not be implemented. Indeed, the following "pure light-blue" move 
clearly satisfies detailed balance: Starting from a spin-configuration s, between all neighbor- 
ing pairs where Si = Sj = H — or Si = Sj = — h independently place a light-blue bond 
with probability 1 — e~ 4/3Ji J . Clusters are identified and, independently, each cluster is flipped 
or left alone with probability \. Such moves (which are obviously not ergodic) can then be 
sandwiched between single-spin updates. 

In addition (for Ising systems) let us recall that it is only the statistical behavior of single 
layers that is actually of interest. Thus we may consider an ensemble of replicas where, in the 
cluster moves, members of the ensemble are paired randomly. Furthermore, if there is some 
translation invariance (e.g. the staggered-field problem) then the replicas may be translated 
relative to one another. Of course neither of the above moves will restore ergodicity; other 
types of updates are still required. 

A closely related set of ideas described in [HBb] and [HBc] has recently been brought to 
our attention. Here the idea is to "fold" the lattice through an axis of symmetry; e.g. the 
xi-axis. Spins at the sites i = (xi,x%, . . .Xd) and i* = (—xi,X2,---Xd) are now regarded 
as "pairs" Assuming that hi = hi* (as would be the case for the staggered field problem) 
one may implement a procedure to generate light-blue clusters as described above. In the 
references cited, single cluster methods were used. Aside from the "boundary spins" , i.e. the 
xi-axis, the method of folding and of duplication are, for all intents and purposes, identical 
for problems with periodically varying fields. However, the advantage of having two separate 
replicas is manifestly apparent for systems where hi is random or, more generally, has no 
spatial periodicity. 

Recently, the authors tested some of these strategies on the square lattice staggered-field 
problem [RMC] . The algorithm was applied to a two- replica system and consisted of light-blue 
cluster moves, random translations of one replica relative to the other by even or odd lattice 
vectors, and Metropolis updates. In the case of odd translations, one replica is also globally 
flipped relative to the other. In [RMC] we use the terminology "active sites" ("inactive sites") 
to refer to sites in light-blue (dark-blue) clusters since, in the cluster part of the algorithm, 
only active sites are flipped. The results were quite satisfactory; for systems up to scale 256, 
our data is consistent with a critical slowing exponent less than 0.5 (compared with local 
dynamics where the value is found to be in excess of 2). 
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